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We study resonant response of an underdamped nanomechanical resonator with fluctuating fre- 
quency. The fluctuations are due to diffusion of molecules or microparticles along the resonator. 
They lead to broadening and change of shape of the oscillator spectrum. The spectrum is found 
for the diffusion confined to a small part of the resonator and where it occurs along the whole 
nanobeam. The analysis is based on extending to the continuous limit, and appropriately modify- 
ing, the method of interfering partial spectra. We establish the conditions of applicability of the 
fluctuation-dissipation relations between the susceptibility and the power spectrum. We also find 
where the effect of frequency fluctuations can be described by a convolution of the spectra without 
these fluctuations and with them as the only source of the spectral broadening. 
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I. INTRODUCTION 

Nano-mechanical resonators are attracting interest in 
various areas of physics. Because they are small and their 
vibrations can be strongly underdamped, even a small 
perturbation can lead to a detectable change of their fre- 
quency. This can be used for charge^^ and mass^j^Jsens- 
ing, high resolution magnetic force microscopyJ^ESl anc j 
other measurements, see Refs. ITT1 - [T51 and papers cited 
therein. The analysis of the frequency change usually re- 
lies on the assumption that the properties of the system 
do not change during the measurement. For example, 
in mass sensing it is assumed that the massive particle 
attached to the resonator does not move. The motion 
would lead to variations of the vibration frequency in 
time. This is because the resonator displacement in the 
vibrational mode depends on coordinates, for example, 
for the fundamental mode of a doubly clamped nanobeam 
it is maximal at the center, whereas for a cantilever it is 
maximal at the apex. The larger the displacement at the 
particle location the stronger is the particle-induced fre- 
quency change.^ If the motion is random, there emerge 
frequency fluctuations, which broaden the spectrum of 
the resonant response. 

The effect of frequency fluctuations on the spectrum of 
an oscillator has been well understood in the limit where 
the correlation time of such fluctuations t c is small. In 
the limit t c — > 0, if the fluctuations are Gaussian, they 
lead to diffusion of the oscillator phase, keeping the oscil- 
lator power spectrum Lorentzian, cf. Ref. [TTJ The effect 
of Gaussian fluctuations with a finite correlation time 
has been also discussed in the literature, see Ref. fl8l - 
and references therein. In the context of nanores- 
onators, an important role can be played also by ran- 
dom freque ncy jum ps due to molecule attachment and 
detachmentSEflEHl 

In the present paper we consider the effect of phase 
fluctuations due to continuous in time random frequency 
variations, which are generally non-Gaussian. A simple 
physical mechanism of such variations is diffusion of a 



massive particle along a nanoresonator. We develop a 
general method for describing the susceptibility of the 
oscillator with continuously fluctuating frequency. One 
might think that this susceptibility could be described 
by weighting a partial susceptibility for a given frequency 
(with the imaginary part described by a Lorentzian) with 
the probability density to have such frequency. How- 
ever, susceptibilities with close frequencies interfere; in 
other words, to find the susceptibility one should add 
the amplitudes rather than the cross-sections of the cor- 
responding transitions. We develop a method that takes 
this interference into account. We then apply the results 
to models of interest for particles diffusing along nano- 
resonators. 

Another question of interest is the interrelation be- 
tween the power spectra and the susceptibilities of under- 
damped oscillators in the presence of nonequilibrium fre- 
quency fluctuations. We provide the conditions of appli- 
cability of the standard fluctuation-dissipation relation, 
including the case of oscillators with weakly nonlinear 
restoring force and nonlinear friction. We also address 
the question of where the effects of oscillator decay and 
thermal fluctuations, on the one side, and of its frequency 
fluctuations, on the other side, can be considered inde- 
pendently. This analysis provides a link to the classical 
work on the lineshape of magnetic resonance in the pres- 
ence of transition frequency fluctuations, where different 
methods were developedPfflSU 

In Sec. |TT] we describe the model, a single mode res- 
onator with a massive particle diffusing along it. We as- 
sume that the vibrations do not affect the diffusion, that 
is, there is no back-action. We introduce the concept of 
partial susceptibility density (PSD) for a given particle 
position and find it in the limiting cases of slow and fast 
frequency fluctuations. In Sec. |III| we derive an equation 
for the PSD of an underdamped oscillator. This equation 



is solved in explicit form in Sec. IV for diffusion confined 
to a small part of the nanoresonator. In Sec. [V] the PSD 
is found in the form of a continued fraction for diffusion 
along a doubly-clamped nanobeam. In Sec. |VI| we study 
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the connection between the susceptibility and the power 
spectrum and find the conditions where the averaging 
over frequency fluctuations can be done separately from 
the averaging over thermal fluctuations of the oscillator. 



We assume that ujq is the largest frequency in the system, 



In Sec. VII we provide a summary of the results. 



II. UNDERDAMPED OSCILLATOR WITH A 
DIFFUSION-MODULATED FREQUENCY 

Mechanical nanoresonators typically have a well- 
separated fundamental mode with eigenfrequency ujq that 
largely exceeds the decay rate T, with the Q factor 
Q = lj /2T ~ 10 3 - 10 5 , se e! 2 ' 3 ' 7 ' 8 ' and references therein. 
Forced vibrations of such mode can be described by the 
model of a driven oscillator, and for not too strong driv- 
ing the oscillator can be assumed harmonic. We can then 
write the equation of motion for the oscillator coordinate 
q in the form 

q + 2Tq+[cj + A D {x)] 2 q = 2— coswt + fr(t). (1) 

Here, F and to are the amplitude and frequency of the 
driving force. The term £t(£) represents thermal noise, 
and M is the oscillator mass. 

The term Ao describes frequency fluctuations. They 
can have different physical origin. In this paper we are 
interested in fluctuations which are continuous in time, 
but are not necessarily Gaussian and have a finite cor- 
relation time. As mentioned above, for concreteness we 
assume that they are caused by a particle absorbed on 
the vibrating nanobeam or trapped in the microchan- 
nel inside the vibrating cantilever and diffusing along the 
nanoresonator. Such diffusion changes the mass distri- 
bution and therefore causes a frequency shift. If x is the 
particle position along the nanoresonator, we can write 



A D (x) = u {m/M)K{x), 



(2) 



where m is the mass of the particle, and lZ(x) can be 
called mass responsivity function; it arises because the 
frequency change depends on the relative amplitude of 
the vibrational mode at the location of the particle^. 

The particle diffusion is described by the Langevin 
equation 



x = -d x U(x)+£ D (t). 



(3) 



Here, U(x) is the trapping potential and is a white 
Gaussian noise, (£d(£i)£d(^2)) = 2D6{t\ — t?), where D 
is the diffusion coefficient; (. . .) indicates ensemble aver- 
aging. The potential U(x) can be created by a droplet 
of a "glue" that confines the attached particle to a small 
region on the nanobeam (a functionalized target area); 
alternatively, we will also consider the case where the 
particle is allowed to freely diffuse along the nanobeam. 

We will be interested in the parameter range where the 
reciprocal correlation time t~ x and the standard devia- 
tion A of the fluctuations of Auix) are comparable, 



T,A,t c ,\6u\ <w, 



SU! 



■ wo- 



(4) 



Conditions Q have been essentially used in writing 
Eq. ([lj where we ignored the effect of the mass change 
due to particle diffusion on the appropriately scaled de- 
cay rate of the resonator and the field amplitude. 

Other types of frequency fluctuations and their effect 
on the osci llator spec trum have been studied in several 
contextsP^ 120 1 23 1 27 ^ We will assume that, even though 
the fluctuations are small on average compared to loq, cf. 
Eq. Q, the interrelation between A and T can be arbi- 
trary. Our results are not limited to Gaussian frequency 
fluctuations. We note that our formulation ignores the 
backaction of the oscillator on diffusion. This backaction 
may lead to nontrivial effects like bistability of forced 
vibrations, which will be studied elsewhere. 



A. Resonant susceptibility in the limiting cases 

We will be interested in the oscillator susceptibility 
X(u), which relates the average value of the coordinate 
to the field, 



(q(t)) = X(u)Fe 



(5) 



we assume that (q) = in the absence of driving. We 
note that the ensemble averaging in Eq. ^ should be 
taken with care in the case of a single oscillator. In the 
experiment, the system is usually assumed to be ergodic. 
However, the ergodicity is established over the correla- 
tion time of frequency fluctuations t c , and for the mea- 
surement time shorter than t c the system may be noner- 
godic. 

The shape of X(u>) near resonance, \u — u>o\ ^ is de- 
termined by the interrelation between the oscillator decay 
rate T, the typical frequency dispersion A, and the cor- 
relation time t c . Frequency fluctuations can significantly 
affect the spectrum for T, t^ 1 < A. 

The susceptibility takes a simple form for compara- 
tively large fluctuational frequency spread, t~ l <C A. In 
this case A gives the typical width of the spectrum. The 
limit t c — > oo corresponds to inhomogeneous broadening, 
where there is no averaging of the eigenfrequency due to 
motion of the particle, there is just a probability for the 
oscillator to have different values of the eigenfrequency. 

To zeroth order in t' 1 , the susceptibility is a superposi- 
tion of what can be called partial susceptibilities x(x;w), 
the susceptibilities for instantaneous fixed positions x. 
More precisely, given the continuous character of the un- 
derlying diffusion, x(ir;u;) should be called the partial 
susceptibility density (PSD). For t c — > oo, Im x(x;w) is 
a Lorentzian centered at frequency ujq + An(x), 



A = ([A D (x) - (A D (x))} 



2\l/2 



Af(a)) = (2M W )- 1 xH; X {u) = I dxx(x;u), (6) 
X (x; u) = iP(x) {T - i[5uj - A D (i)]}^ (t c -> oo). 
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Here, P{x) is the probability density for the diffusing 
particle to be at point x. The overall spectrum Im x( w ) 
in this limit is typically non-Lorentzian. It becomes par- 
ticularly simple in the limit r — > 0, in which case 

Im x (w) = irJ2 Xa P(x u )/\A' D (x u )\ 

with x u given by equation Ar>(x u> ) — Sui; A' D (x) = 
d x A D (x). 

In the opposite limit, t^ 1 ^> A, the oscillator can- 
not "resolve" frequency variations, they are averaged 
out. This is similar to the motional narrowing effect in 
NMR. To zeroth order in t c A we expect Im x( w ) to be a 
Lorentzian curve centered at frequency luq + (Ad) with 
halfwidth T, where 



(A D ) = J dxP{x)A D {a 



Clearly, the shape of x( w ) i s qualitatively different in 
the opposite limits of t c A. In what follows we will de- 
velop an approach that allows one to find the suscepti- 
bility for an arbitrary t c A. We will also relate the results 
to the analysis of dcphasing developed by AndersorP^l 
and Kubo and TomitcP^ES] m the context of resonant ab- 
sorption by two-level systems. For diffusion described by 
Eq. Q , the system has detailed balance in the absence of 
periodic driving (again, we disregard the effect of back- 
action on the diffusion), and therefore the susceptibility 
can be obtained from t he p ower spectrum of the oscilla- 
tor calculated for F = 0p2 However, we will calculate the 
susceptibility directly, since our approach applies also to 
systems without detailed balance. 



III. EQUATION FOR THE PARTIAL 
SUSCEPTIBILITY DENSITY 



It is convenient to analyze resonant response of the 
oscillator using the standard method of averaging. This is 
done by changing from q(t), q to slow variables u(t), u*(t), 



q(t) = ue iu>t + u*e- iut , 

'(7) 

The Langevin equations of motion for slow variable u(t) 
in the rotating wave approximation (RWA) is 

u = -[r + i(8w - A D (x))]u + — — + £r ;tt (t) (8) 

ZJMU 

where £t;u(£) = (2iw) _1 £r(i) exp(— itot) is the random 
force. Equation for u* can be obtained from Eq. ([8| 
by complex conjugation. We note that, in fact, the 
Markovian equations of motion for u, u* have a much 
broader range of applicability than the original equation 
([TJ. They apply even where relaxation of the oscillator 
is not described by a simple viscous friction force, as in 
Eq. 0, but is delayed. Quite generally, the delay disap- 
pears on the slow time scale ~ 1/r. The random forces 
£,T-.u{t) and Ct-«(^) are a l so ^ correlated on the slow time 



q = iuj (ue wt - u*e- iut ) 



scale rather than in the "fast" time, see Refs. I3H1 and |3"T1 
and papers cited therein. 

The probability distribution of the oscillator in slow 
time p^UjU* ,x;t) is described by the Fokker-Planck 
equation,^ which follows from Eqs. (|3|) and (|S|) 



d t p = d u ([r + i(5ui - A D (x))} up) - d u p - 

2iMu> 



d u * ([T-i(8u}-A D (x))]u*p) + 



F 



2iMu 



d U 'P- 



Tk B T 



dl u -P + L D [p} 



(9) 



Here, T is the bath temperature and is the diffusion 
operator, 



L D [p]=d x {pd x U) + Dd 2 x p. 



(10) 



The scaled susceptibility xi 1 ^) f° r u close to loq is given 
by the expectation value (2Mu)/F)(u*). It is convenient 
to write it in the form of an integral over x of the PSD 
X(x; u), see Eq. Using Eqs. ([5| and ^ one can write 
the PSD in the form 



X{x;u) 



2Mw 
F 



dudu* u* p s t{u,u* ,x), (If) 



where p st (u,u* ,x) is the stationary solution of Eq. |9|. 
Multiplying Eq. ^ by u* and integrating over u and u* 
one obtains 

[r - i(Su - A D (x))} x(x;oj) - L D [x(x;u))] = iP(x), 
P(x) = f dudu* Pst {u,u*,x) = Z - l e~ ui - x) l D (12) 



with Z = Jdx exp[-U(x)/D}. 

Equations ([6| and ( 12 ) reduce the problem of the spec- 
trum of the oscillator to solving a diffusion equation for 
the PSD x(x;cj). They show also that the values of the 
PSD for different particle positions x are coupled to each 
other. This coupling becomes small for large A, i.e., for 
the case where the actual range of Ajj(x) in Eq. 
is large. If the drift and diffusion operator Ljj[x 
Eq. (12) is disregarded, one immediately obtains Eq 



(12) 



m 
G) 



for x(x; ui). The (minus) reciprocal correlation time t 
is given by the lowest eigenvalue of Ljj, and disregarding 
Ljj is justified provided t c A 3> I. 

In the opposite limit where the variation of Ad (x) with 
x can be disregarded in Eq. ( 12 ), 



X (x;u } )^iP(x)/[T-i(6uj-(A D ))}. 

Formally, one can think that the PSDs have the same 
shape in this case, but it is more correct to say that 
they are strongly coupled by the operator and the 
contributions to x( w ) °f the PSDs from different ranges 
of x cannot be identified. 



The solution of Eq. ( 12 ) can be written in the form 



of a convolution of the "complex Lorentzian" suscepti- 
bility 7r _1 (r — iSuj)^ 1 and the PSD Xph^;^) calculated 
in the absence of the oscillator decay and determined 
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by phase fluctuations only. This form applies also in a 
more complicated case of a nonlinear oscillator as long as 
Ac is independent of the oscillator amplitude, see Sec- 
tion [VTj For a nonlinear oscillator, the additive thermal 
noise £t(*) affects the susceptibility, in contrast to the 
case of a linear oscillator. 

In the general case of an arbitrary t c Ajj, an arbitrary 
form of the frequency shift Ad(x), and an arbitrary po- 
tential U(x), Eq. (12 1 can be solved numerically. How- 



ever, there are important situations where an analytical 
solution can be obtained. They are discussed in the next 
Section. 



IV. HARMONIC CONFINING POTENTIAL 

The susceptibility x(w) can be found in an explicit 
form if the diffusing particle is confined to a small region 
of the nanobeam. A simple and important form of the 
confinement is described by a parabolic potential U (x) — 
k(x — xo) 2 /2. Such potential models confinement due to 
a small droplet of a polymer "glue" on a nanobeam or 
to centrifugal forces created in a suspended nanochannel 
by additional drivingPWe will assume that the potential 
minimum xq is far from the ends of the resonator and that 
the typical displacement of the particle from the potential 
minimum (D/k) 1 ^ 2 <C L. Then the shift Ad(x) can be 
expanded about the value Ad(xq) which we will set equal 
to zero (it can be incorporated into uj ), 



A D (x) w a(x — x ) + (3(x — x ) z 



(13) 



For small x — xq the linear in x — xq term in Eq. ( 13 ) 



generally dominates and the quadratic term can be dis- 
regarded. However, if the linear term is small, it is nec- 
essary to keep the quadratic term. This is the case, 
for example, if Xq is at the center of a double-clamped 
nanobeam. For diffusion in a parabolic potential, the 
correlation time of x(t), and thus of A^, is t c = l/k. 
We seek the susceptibility as 



X{x\(J) = 



(It 



At 8lo ~ 



(14) 



with initial conditions 



A(0) = -k/2D, B(0) = 0, C(0) = - ln(fc/27r£)). 



Equations ( 17) can be easily solved, and the solution is 



expressed in elementary functions. It allows finding the 
function \(t) = J dxx(x; t), 



m = i 

xH = 



A(t) 



1/2 



exp 



B 2 (t) 
4A(t) 



C(t) 



dte itSw x{t). 



(18) 



Since the general solution is somewhat cumbersome, we 
will consider the cases where only one of the coefficients 
a and (3 is nonzero. 



A. Frequency change linear in the particle 
displacement 

We start with the case Ad(x) = a(x — xo). For U(x) — 
kx 2 /2 this case corresponds to the oscillator frequency 
being modulated by the Ornstcin-Uhlenbeck noise, which 
is an exponentially correlated Gaussian noise. Indeed, 
from Eq. ([3| x is Gaussian, and thus Ad(x) cx x — xa is 
Gaussian, too, with (A^) = and 



(A D (x(t))A D (x(t'))) = AV'!"'! 



with A = {a 2 DJkf 12 
m E( 

X(t) = iexp 



From Eqs. (|17|) and (j 1 S|) for j3 — we obtain 

-kt\ 



A 2 A 2 , 



(19) 



(20) 



As explained in Section |VT] below, this result could be 
also obtained from the expression for the power spectrum 
of the oscillator without driving by using the cumulant 
expansion. Equation ( 20 ) is equivalent to the result of 
Ref . M 



It is interesting to compare Eq. ( 20 1 with the asymp- 
totic results for small and large t c A A/k. For t c A ^ 1, 



one can expand the exponent in Eq. (20) in kt, getting 



From Eqs. (|12| and (|14|) we obtain an equation for func- 

L D [ X }=0 (15) 



X(x; 0) = iP(x) = i(k/2irD) 1 / 2 exp[-fc(a; - x Q ) 2 /2D}. *(*) « ex Ptr R " A ^ 2 /2)- Substituting this expression 

into Eq. ( 18 ) one obtains an expression that coincides 
with Eq. with the corresponding P(x) and Ar>(x). 
Specifically, Im x(<j-0 is a convolution of a Lorentzian with 
halfwidth T and a T-independent Gaussian distribution 
cx exp(— <5w 2 /2A 2 ). In the opposite limit, t c A <C 1, the 
peak of Im xi 10 ) IS Lorentzian, with halfwidth Y + t c A 2 . 
As a whole, the spectrum of Im x( w ) is symmetric, with 
maximum at w =wo- 



tion xi x 'i t) of the form 

d tX +[r + iA D (x)]x- 



with x(^;0) given by Eq. (14). 

Equations (fT3|) and |l5|) have a solution 



X (x; t) = i exp [A(t)(x - x ) 2 + B(t)(x - x Q ) + C(t)) 

,( 16 ) 

where functions A, B, C are given by a set of ordinary 
differential equations 



.4 

C 



ADA 2 
D(B 2 



2k A 
2A) - 



-0, 
k-T 



B = kB + ADAB 



(17) 



B. Frequency change quadratic in the particle 
displacement 

We now consider the case where A D {x) — f3(x — 
xq) 2 . As mentioned before, this case is interesting if 
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the equilibrium position of the particle is at an antin- 
ode of the vibrational mode of the nanomechanical res- 
onator. The value of (3 is easy to estimate using the 
standard analysis.^ For example, for the fundamental 
mode of a doubly clamped beam of length L we have 
P = ujomir 2 1 AL 2 M . Frequency fluctuations due to dif- 
fusion are non-Gaussian, with (A^) = (3D/k and A = 
y/(A 2 D ) - (A D ) 2 = V2\(3\D/k. The correlation time is 
the same as for the linear in x frequency change, t c = l/k. 
From Eqs. ((171) and (Il8j) we find 



X(t) = 2iaJ 2 exp 



-Tt - -k(ap 



Dt 



x [(a fj + l) 2 - (ap - l) 2 exp(-2apkt)] 
ap = (1 + AipD/k 2 ) 1 ' 2 (Re ap > 0) 



-1/2 



(21) 



This expression shows that, for the frequency shift 
quadratic in the particle displacement, decay of \(t) is 
non-exponential in time, which means that the spectrum 
of Im x(cj) is non-Lorentzian. As expected from the gen- 
eral arguments, except for the trivial factor exp(— Yt) 
that describes decay of the vibration amplitude, x(t) is a 
function of the scaled time kt = t/t c and one dimension- 
less parameter ap = (1 + Ai(Ao) /k) 1 ^ 2 ; this parameter, 
in turn is given by the ratio of the standard deviation of 
the fluctuating frequency A = \/2|(A£>)| to the recipro- 
cal correlation time of fluctuations k. 



From Eqs. ( 18 1 and ( 21 1 , for A <C k, the major effect of 
diffusion is the shift of the peak of Im x( w ) by « (Ad); 
to the first order in A/k the peak of Im x( w ) remains 
symmetric and Lorentzian. However, for large A/k the 
peak becomes strongly asymmetric. For arbitrary A/k 
one can write x( w ) as 

1/2 ^ (2n-l)H (aff-1) 2 " 

XM=2a/ g 2 n n[ + 1)2^1 *nM, 



Xn{u) = t{T+-k[{An+l)a p -l}-i5uj} . (22) 



Equation (22) presents the susceptibility in the form of 
a sum of the spectra % n (o;). Functions Im x„(cj) have 
peaks at equally spaced frequencies ojo+(4n+l)fclm ap /2, 
with halfwidth V + [(An + l)Re ap - l]k/2 that lin- 
early increases with n. We note that these spectra 
should not be called partial spectra of the oscillator; even 
in the limit A 3> k the distance between their peaks 

2 5 / 4 (£;A) 1 / 2 is generally smaller than their halfwidth 
~ 2~ 3/4 (4n + 1)(A:A) 1/2 + Y. Moreover, functions x„(w) 
enter the expression for x( w ) with complex weighting fac- 
tors, so that Im x( w ) is determined by the both real and 
imaginary parts of ^n(w). 

Equation ( |22[ ) is convenient for a numerical evaluation 
of x( tJ )- K also allows establishing the connection with 
the limit t c A = A/k > 1, Eq. Q. To do this one no- 
tices that \ap\ 3> 1 for A ^> k. Typical values of 8u> 
within the peak of Im x( w ) are ~ @D/k = (Ap) [see 
Eq. (pi below]. Since [P[D/k 2 w \ap\ 2 /A > \ap\, the 




8co/A 



FIG. 1. The scaled susceptibility Im x( u ) f° r the case of 
dephasing due to a particle diffusing about the antinode of 
a nanoresonator, with the frequency shift quadratic in the 
particle displacement. The data refer to the resonator energy 
decay rate r = 0. The frequency 8uj — uj — wo is scaled 
by the standard deviation of the oscillator frequency due to 
the massive particle A. Curves 1 through 5 refer to t c A = 
100, 40, 3, 0.05, and 0.005 respectively. 



major contribution to the sum over n in Eq. (22) comes 



from n 3> 1 and one can replace summation over n by 
integration. The integrand has a singularity for n — n p , 
where n p = iSuj/2kap for T = 0, \n p \ 3> 1 for typical 
Slu. Integration over n can be done by lifting or lowering, 
depending on the sign of /3, the integration contour up 
to Im rip, which gives for T = 

Im xM pa (7rk/2pDSuj)^ 2 e- kSu]/2pD e{l3Sio), (23) 

where 0(x) is the step function. Equation ( p3| applies 
for 



A»fc,r, 



|<M » (kA) 1 ' 2 . 



The spectrum ( 23 ) has a very specific shape that makes 
it possible to identify the corresponding mechanism 



It 



is profoundly asymmetric, with a square-root divergence 
near the maximum in the neglect of corrections cx k/A 
and with an exponential tail. 

The evolution of the spectrum Im x( w ) with varying 
A/k is seen in Fig. [I] For convenience, the figure is 
plotted for Y = 0; to allow for Y the spectra should be 
convoluted with the Lorentzian distribution. The sus- 
ceptibility as a function of dimcnsionless frequency Suj/k 
depends on the single dimensionless parameter A/k. It 
is seen from Fig. [l] that, with increasing A/k, the spec- 
trum shape changes from an almost symmetric peak cen- 
tered close to (Ad) with width cx A 2 /fc for small A/k 
to the strongly asymmetric distribution that approaches 
Eq. pi for large A/k. 



V. UNCONFINED DIFFUSION ALONG A 
NANOMECHANICAL RESONATOR 

We now consider dephasing due to a particle that freely 
diffuses along a nanoresonator, but does not leave it. We 
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assume that the nanoresonator is a one-dimensional dou- 
bly clamped beam, and we are interested in its funda- 
mental mode. The change of the mode frequency due 
to a particle at a point x is determined by the squared 
vibration amplitude at x.^ For a beam of length L this 
gives An(x) = — 7 cos 2 (irx/L), where 7 = luqiti/M: x is 
counted off from the center of the beam. The stationary 
probability distribution of the particle along the beam is 
uniform, P(x) = l/L. 

The average frequency shift and the standard fre- 
quency deviation are, respectively, (An) — —7/2 and 
A = j/y/8. The correlation time of frequency fluctu- 
ations can be found by calculating the time correlation 
function of An(x(t)), which can be done following the 
standard prescription^! (the classical analog of the quan- 
tum regression theorem) . It involves evaluating the prob- 
ability density pa of a transition (xo,t — 0) — > (x, t) 
integrated over xq with the appropriate weight. Func- 
tion pa is given by the solution of diffusion equation 
PA(x.t) = —Dd%pA(x;t). The boundary conditions fol- 
low from the absence of current, d x pA — for x — ±L/2, 
and the initial condition for the correlation function of 
A D is p A (x;t = 0) = [A D (x) - (A D )]/L. This gives 

(A D (x{t))[A D (x{Q)) ~ (A D )]) = [ dxA D (x)p A (x;t) 



A 2 exp(-t/t c ), 



t: 



D(2ir/L) 2 



(24) 



The oscillator susceptibility is given by Eq. ( 12 ) with 
U(x) — and with boundary conditions d x x(x;uj) = 



for x = ±L/2. It is clear from the structure of Eq. (12) 
and the expression for An(x) that the solution can be 
sought in the form 

x( x ] w ) = b n (ui) cos(2i:nx / L) . 

n>0 



Equation ( 12 ) is then reduced to a tri-diagonal system of 



linear equations for coefficients b n , which can be solved 
by the method of continued fractions. This gives for x(w) 



^0 



A 2 



(25) 



Vi 



A 2 /2 



V 2 



A 2 /2 



where V n = T + nH' 1 - i(6u - (A D )). We note that, 
alternatively, x{ u ) can be also expressed in terms of the 
Mathieu functions. 



From Eq. (251, for V = the reduced susceptibility 



x{oj)/t c is a function of dimensionless frequency t c 8uj that 
depends on the single parameter, t c A. For t c A <C 1, to 
find the peak of Im x( w ) one can ignore in Eq. (25 1 frac- 
tions that contain V n with n > 2. This gives a Lorentzian 
peak, 

Im xM « A 2 t c /[{8u - (A D )) 2 + (A 2 t c ) 2 ] 




5to/A 



FIG. 2. The scaled susceptibility Im x^) f° r the fundamen- 
tal mode of a doubly clamped resonator with a particle freely 
diffusing along it. The data refer to the resonator energy de- 
cay rate F = 0. The frequency 8ui — ui — ujq is scaled by the 
standard frequency deviation A. Curves 1 through 5 refer to 
t c A = 50, 5, 0.5, 0.05 and 0.005, respectively. 



for T = 0. The halfwidth of the peak A 2 t c has a typi- 
cal form of the width of the spectral peak for motional 
narrowing in NMR, see Sec. |VI| 

In the opposite limit, t c A ^> 1, we obtain from Eq. ^ 
or Eq. (p5l for T = 



Im xM « 2 [< Ad > 



(8lj — (An)) 2 



-1/2 



This spectrum as a function of u> has two inverse square- 
root peaks symmetrically spaced around frequency u>o + 
(Ad). 



Expression (251 is convenient for a numerical evalua- 



tion of the susceptibility in the general case of arbitrary 
t c A. The evolution of the shape of Im x( w ) with varying 
t c A for r = is illustrated in Fig. [2] The spectrum re- 
mains symmetric, but as expected from the asymptotic 
expressions, can have a single peak or two peaks. 



VI. FLUCTUATION-DISSIPATION RELATION 
AND THE SEPARATION OF PHASE 
AVERAGING 

Along with the susceptibility, the power spectra of 
nano- or micromechanical resonators are also often stud- 
ied in the experiment!^ For systems in thermal equi- 
librium, the two spectra are simply related by the 
fluctuation-dissipation theorem. However, frequency 
modulation by an attached diffusing particle (or by an- 
other external source) can drive the system away from 
equilibrium. Here we derive the conditions where the 
fluctuation-dissipation relations apply in the presence of 
frequency modulation. Another important issue that we 
address is whether it is possible to perform averaging over 
thermal fluctuations and over the externally imposed fre- 
quency fluctuations independently. 

We will be interested in an underdamped oscillator. It 
is characterized by slow variables u,u* , Eq. ([7|. Fluctu- 
ations of these variables in slow time are usually Marko- 
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vian, in the RWA, and can be described by the Fokker- 
Planck equation, cf. Eq. ([9|. If the restoring force of 
the oscillator is weakly nonlinear, i.e., the oscillator po- 
tential is weakly nonparabolic, the major effect of this 
nonlinearity is that the oscillator frequency depends on 
the vibration amplitude. Also, if the friction force is non- 
linear, the friction coefficient is amplitude-dependent; in 
this latter case the form of the operator that describes 
the effect of the thermal noise associated with friction 
changes compared to Eq. Q. 

In the absence of diffusion of an attached particle the 
Fokker-Planck equation for a weakly nonlinear oscillator 
in slow variables was derived e arlier for both microscopic 
and phenomenological models Because of the diffu- 
sion the oscillator parameters slowly vary in time. On 
the other hand, the diffusion itself may depend on the 
vibration amplitude. 

For the analysis below it is convenient to introduce the 
slow variables in such a way that they are independent 
of the drive frequency, 

q(t) = uq exp(iwoi) + Uq exp(— iwot), 

q = iujQ [uq exp(iwoi) — Uq exp(— iojQt)] . (26) 

In these variables the Fokker-Planck equation reads 

dtp = L FP [p] - [{F/2iMu:)e lt ^d Uo p + c.c] , (27) 

where 



L F p[p] = [d Uo (K(\uo\ 2 ,x)u p) 
, k B T 



c.c. 



d 2 Uo< [T nl (\u \ 2 ,x)p\+L D [p] (28) 



In the RWA, functions K and r„; depend only on the 
scaled squared vibration amplitude |ito| 2 , but not on uq 
and Uq taken separately! 30 * 31 ^ This can be understood 



from Eq. (26): prior to averaging over the period 2tt/ujq 
in the RWA, the corresponding terms would be propor- 
tional to fast oscillating factors exp(±iu;o£), and therefore 
in the RWA they average to zero. The real part of K 
describes dissipation due to coupling to a thermal reser- 
voir, whereas the term oc k B TT n i describes fluctuations 
induced by the reservoir. 

Functions K and r n j can also parametrically depend 
on the particle position x. We assume that this depen- 
dence is such that the fluctuation-dissipation relation be- 
tween K and r„; holds for any x, 



Re K(r,x) - T nl (r,x) + 



k B T 
2Muj 2 



drTnl = 0, 



(29) 



where 



The diffusion operator Lp, can depend on the \uq\ . We 



note that in Eq. (27), in contrast to the Fokker-Planck 
equation ([9|, K does not depend on u>, but the term oc F 
has a time-dependent factor. Respectively, the solutions 
of Eqs. ^ and Q for p are also different even if we use 
the same model of the oscillator as in Eq. ([9]); however, 
they are connected by a simple canonical transformation. 



In the absence of modulation the po wer s pectrum of 
the oscillator near resonance (u> « ujq) i d 31 * 32 ! 

/>oo 

Q(u) = Tr^Re / dte iut {q(t)q(0)) F = 
Jo 

WTT^Rey dte ttS " J du du* dxu* oPuo (t). (30) 

Here, Suj = u - uj , p Uo (t) = p U0 {u Q ,u* Q ,x;t). Func- 
tion p Uo (t) is given by Eq. |27|) with F = 0. The 



initial condition is p Uo (t 



uopeq, where p eq 



Pcq(wo, Uq, x) is the equilibrium distribution for F = 0. 
Formally, from Eq. (l27|) one can write 



L F p[p cq ] = 0, (31) 



Pu {t) = exp(L F pt)u p e 



with operator L F p given by Eq. (28). This, combined 



with Eq. ( 30 1 , provides a formal expression for the oscil- 



lator power spectrum. 

On the other hand, from Eqs. §5y, and (27), res- 
onant susceptibility xi 1 ^) is determined by the average 
value of Uo(i) which, to first order in F, is given by the 
linearized solution of Eq. (27 1, 



dte liSu) I du Q du*Qdxu* 



x exp(L FP t)d U 'p cq . 



(32) 



The fluctuation dissipation relation for the scaled sus- 
ceptibility means that, near resonance, there should hold 

Im xH = (2nMuj^/k B T)Re Q(oj). 



A comparison of Eqs. (30) and (31 ), on the one side, and 



Eq. (32 1, on the other side, shows that this relation ap- 
plies if d u *p cq = -C T u p cq , with C T = 2Mu>l/k B T. In 
turn, this condition holds if functions Re K and r„; are 



related by the fluctuation-dissipation theorem Eq. (29) 
and operator L B is independent of |uo| 2 - Indeed, from 



Eq. (29) it follows that, for fixed x, the equilibrium dis- 



tribution of the oscillator over uq, Uq is of the Boltzmann 
form, p cq cx p b (u ,Uq), 



Pb(u ,u*q) 



Mujq 1 



cM-2Mu; 2 Q\uo\ 2 /k B T). (33) 



Here, 2Mu} 2 \u \ 2 = | (A/wgg 2 + Mq 2 ) is just the oscilla- 
tor energy, neglecting small nonlinear corrections. When 
Lb is independent of |wo| 2 , the equilibrium distribution 
over x is determined by a factor P(x), L B [P] — 0; for 
the model of diffusion used in this paper P(x) is given 
by Eq. ( 12 ). The whole equilibrium distribution is multi- 



plicative, it is a product of functions P(x) and p b {uq, Uq) 
that depend on x and uq,Uq separately . 



A. Convolution representation 

The calculation is significantly simplified in the im- 
portant case where the susceptibility can be written as a 



convolution, 



From Eqs. (34), (35), (37), and the fluctuation- 



XH = / du'x oac {uj')xD{u - u>'), 



Xd{u) = / dxxD{x;uj). 



(34) 



Here, Xosc(^) is the susceptibility in the absence of diffu- 
sion and xd(x;co) is the susceptibility that describes the 
effect of diffusion independently from the oscillator dy- 
namics. This representation applies, in particular, for the 
models of frequency fluctuations discussed in this paper, 
with x osc (w) = (oj — ljq — iT)^ 1 and with xd calculated 
from Eq. (pi) for T = 0. 



The representation (34) is particularly helpful if there 



holds the fluctuation-dissipation relation between x(cj) 
and Q(w), which allows using the power spectrum to 
find the susceptibility. Still, the very applicability of the 
fluctuation-dissipation relations does not guarantee that 
Eq. ( 34 ) would apply. We now provide the sufficient con- 



dition. 

Finding Q(uj) requires solving the Fokker-Planck equa- 



tion dtPu = Lpp[p Uo ]. From Eq. (28) one can see that 



the solution can be sought in the form 

Pu (u ,u* ,x;t) = u p{\u Q \ 2 ,x;t). 



Equation (34) will apply if function p is a product, 

p(\u \ 2 ,x;t) = p osc (\u Q \ 2 ;t)p D (x;t), (35) 

i.e., in the equation that follows from the Fokker-Planck 
equation for p one can separate variables |it | 2 and x. A 
straightforward analysis shows that this happens if 



d x T nl = Re d x K = 0, Im d 2 x K(r, x) = 0. 



(36) 



In other terms, r„; and Re K should be independent of x, 
whereas Im K should be a sum of terms that depend on 
x and |uo| 2 separately, Im K = Im K osc (\uo\ 2 ) — Ad(i). 
These conditions hold in the model discussed in the main 
part of the paper. We note that, for a nonlinear oscillator, 

Xosc 

(oj) is non-Lorentzian and can be asymmetric!^ 



B. Relation to dephasing in two-level systems 

In the case where the oscillator dynamical variables 
separate from the coordinate of the diffusing parti- 
cle, function pn(x;t) is given by equation dtpD = 
—iAxj(x)pD+Li)\pr)\- A formal solution of this equation 
is 



p D (x;t) = j dxiP(xi)p D (x;t\xi;0) (37) 
p D (x;t\x(Q);0) = /e -i # dt ' A ° ( x(t,) ) 5(x(t) - x) 



where x(t) is given by the Langevin equation ^ and the 
averaging is done over realizations of the noise £d (t) that 
drives the diffusing particle. 



dissipation relation we obtain 



die 1 



itSui I -i f* dt'An (»(«')) 



(38) 



where the averaging is now done both over the realiza- 
tions of£o(t) and over the stationary distribution of x(0). 
Equation ( 38 ) has the same form as the expression for 



the susceptibility of a two-level system with fluctuating 
frequency that was studied in th e cele brated papers by 
AndersorPS and Kubo and Tomita^^SI assuming that the 
system was in thermal equilibrium. In particular, the av- 
eraging in Eq. ( 38 ) is simplified if frequency fluctuations 
are Gaussian, in which case one can use the cumulant ex- 
pansion. This is the case for diffusion in a parabolic po- 
tential with dAp,/dx = const, where the frequency fluc- 
tuations correspond to the Ornstein-Uhlenbeck noised. 
The methods of Refs. EH EMID (see also Refs.|34]) do not 
immediately apply to other cases studied in this paper. 
As we showed in Sees. II-V, in all cases of interest the 
solution is naturally obtained using the method of cou- 
pled partial susceptibilities. We note that this method 
applies also if the system is far from thermal equilibrium 
and the fluctuation-dissipation relation between the sus- 
ceptibility and the power spectrum does not hold. 



VII. CONCLUSIONS 

We have studied resonant susceptibility of an un- 
derdamped oscillator whose eigenfrequency continuously 
fluctuates in time. The analysis is based on the method 
of partial spectral density. Such density corresponds to 
a given eigenfrequency value in the limit of very slow 
eigenfrequency variations. The variations lead, on the 
one hand, to the finite lifetime of states with different 
eigenfrequencies, and, on the other hand, to the interfer- 
ence of the spectral densities for close eigenfrequencies. 
The resulting overall spectrum depends on the interre- 
lation between the bandwidth A of the eigenfrequency 
variations and the correlation time of these variations t c . 

We have developed a method that allowed us to study 
the susceptibility for an arbitrary t c A. It involves de- 
riving and solving a differential equation for the partial 
spectral density. The specific results are formulated for 
nano-mechanical resonators whose frequency can fluctu- 
ate if they have particles diffusing along them and thus 
changing their effective mass. 

Explicit results have been obtained for three models: 
(i) a particle diffusing in a small region centered at a 
general position on the nanoresonator; (ii) a particle dif- 
fusing about the antinode of the vibrational node, and 
(iii) a particle uniformly diffusing along the nanobeam. 
The shape of the absorption peak Im x( w ) is different in 
all these cases, varying from symmetric non-Lorentzian 
in (i), to asymmetric in (ii), to symmetric but possi- 
bly double-peaked in (iii). In all these cases the shape 
strongly depends on the interrelation between A and t~ . 
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Another general result refers to the interrelation be- 
tween the oscillator susceptibility and the power spec- 
trum. We have found the conditions where the standard 
fluctuation-dissipation relation applies in the presence 
of phase fluctuations even where these fluctuations are 
nonequilibrium. In addition, we have established where 
the spectrum of a generally nonlinear underdamped os- 
cillator is a convolution of the spectrum in the absence 
of phase fluctuations and the spectrum broadened by 
phase fluctuations only. The latter broadening can be 
also described, at least in principle, using the meth ods 
developed by AndersorPS and Kubo and Tomita^SEni for 
two-level systems with a fluctuating transition frequency. 
These methods are particularly convenient where the fre- 
quency fluctuations are Gaussian, and our results for 
the case (i) above are equivalent to those obtained us- 
ing thempS The fluctuations discussed in the cases (ii) 
and (iii) are non-Gaussian and have not been previously 
studied, to the best of our knowledge, nor the method of 
coupled partial spectral densities has been used. 

The results of the paper have immediate relation to 
mass sensing with nanoresonators. For the particle that 
is being analyzed and that diffuses along a nanoresonator, 
parameter A is proportional to the particle mass, whereas 
t c is determined by either the particle confinement, as in 
cases (i) and (ii) above, or is inversely proportional to the 
diffusion coefficient D, as in case (iii). The shape of the 
spectrum provides an important additional information 
about the attached particle and its dynamics, compared 
to the conventionally considered shift of the spectral line. 

Observation of the effects of particle dynamics is possi- 
ble for comparatively large diffusion coefficients. Fast dif- 
fusion can happen along carbon nanotubes^ or for par- 



ticles inside low-viscosity nanochannels embedded into 
cantilevers^^ or on solid-state nanobeams at elevated 
temperatures. In particular, for carbon- nanotube based 
nanoresonators of length ~ I /iirPS we get from Eq. ( 24 ) 



the correlation time t c < 10~ 5 s for D ~ 10~ 4 cm ^s 
Such D is smaller than the calculated values of the 
diffusion coefficients for different simple molecules, see 
Refs. and [37J This suggests using spectral measure- 
ments of nanoresonators to determine the diffusion coef- 
ficient in carbon nanotubes. This makes it also possible 
to use temperature as an additional means of the analysis 
of diffusion in nanoresonators. 

After this paper was completed, we learn of the work 
by Yang et alpSI where phase fluctuations due to diffusion 
of particles along a nanoresonator were observed already 
for T < 80 K (in contrast to the present paper, the dif- 
fusion was not confined to the nanoresonator itself, and 
there was an influx of particles to keep their mean num- 
ber constant) . We note that fluctuations of the nanores- 
onator frequency can be due to other reasons, for exam- 
ple, to fluctuations of the charge on the substrate above 
which the nanoresonator i s lo cated or to charge fluctua- 
tions in the nanoresonator P^l The analysis of this paper 
can be extended to this case. 
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